Kondo effect of an adatom in graphene and its scanning tunneling spectroscopy 
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We study the Kondo effect of a single magnetic adatom on the surface of graphene. It was shown 
that the unique linear dispersion relation near the Dirac points in graphene makes it more easy to 
form the local magnetic moment, which simply means that the Kondo resonance can be observed 
in a more wider parameter region than in the metallic host. The result indicates that the Kondo 
resonance indeed can form ranged from the Kondo regime, to the mixed valence, even to the empty 
orbital regime. While the Kondo resonance displays as a sharp peak in the first regime, it has a 
peak-dip structure and/or an anti-resonance in the remaining two regimes, which result from the 
Fano resonance due to the significant background leaded by dramatically broadening of the impurity 
level in graphene. We also study the scanning tunneling microscopy (STM) spectra of the adatom 
and they show obvious particle-hole asymmetry when the chemical potential is tuned by the gate 
voltages applied to the graphene. Finally, we explore the influence of the direct tunneling channel 
between the STM tip and the graphene on the Kondo resonance and find that the lineshape of 
the Kondo resonance is unaffected, which can be attributed to unusual large asymmetry factor in 
graphene. Our study indicates that the graphene is an ideal platform to study systematically the 
Kondo physics and these results are useful to further stimulate the relevant experimental studies on 
the system. 

PACS numbers: 72.10.Fk, 73.20.Hb, 73.23.-b 



I. INTRODUCTION 

The experimental realization of grapheneiS composed 
of a monolayer of carbon atoms trigged a new wave to 
study the carbon-based materials, both from fundamen- 
tal physics and applications^— The graphene possesses 
perfect two-dimensional massless Dirac fermion behav- 
iors, and the valence and conduction bands touch at two 
inequivalent Dirac points K_ and K + at the corner of the 
Brillouin zone. Around the Dirac points the low-energy 
excitations are linear and the unique electronic structure 
leads to a number of unusual electronic correlation and 
transport properties of graphene^r— 

For example, the localized magnetic moment on 
adatom with inner shell electrons in graphene behaves 
quite different to that in metals, a conventional Fermi liq- 
uid with a constant density of states near the Fermi sur- 
face. Anderson showed that in such a host the strongly 
interacting impurity ion can be in a magnetic state if 
the on-site interaction U is sufficiently large and/or the 
hybridization V between the impurity ion and the con- 
duction electrons is sufficiently small^ Explicitly, it is 
required that the energy of the single occupancy states 
Sd is below the Fermi level /i but the energy of the dou- 
ble occupancy states, namely, Ed + U, should be larger 
than fj,. Moreover, the magnetic regime was found to be 
symmetric around y — (/x — Ed)/U = 0.5. For details one 
can refer to the phase diagram given by the well-known 
Ref. [12]. In contrast, when the host is graphene, the 
situation is dramatically different due to the linear den- 
sity of states near the Dirac points Firstly, the mag- 
netic phase diagram is not symmetric, it depends on that 



the impurity level e<j is above or below the Dirac points. 
Secondly and also most important difference is that the 
phase boundary is not symmetric around y — 0.5 in both 
cases. When Ed is above the Dirac points, the magnetic 
state extends to y < 0, which means that even the impu- 
rity level is located above the Fermi energy, it is possible 
that the localized magnetic moment exists. This is quite 
different to that in usual metallic host. In the other case 
that Ed is below the Dirac points, the phase boundary 
extends to y > 1, which means that even the impurity 
is in the double occupancy state, it can be magnetized. 
It was realized that all these unusual features are due to 
the linear dispersion of the graphene around the Dirac 
points, and the impurity level is broadened dramatically 
by hybridization to cross the Fermi level^ 

A direct consequence of the much wide magnetic 
regime is that the Kondo effect in graphene^ may ex- 
ist in a more wide parameter regime than that in the 
usual metallic host in which the Kondo effect only exists 
in the Kondo regime, namely, e^/r <C —0.5^ where T 
is the effective tunneling coupling between the impurity 
and conduction electrons. Moreover, it is also interest- 
ing to explore what effect the dramatic broadening of the 
impurity level on the Kondo physics. This motivates us 
to study the Kondo effect of an adatom on the surface of 
the graphene sheet in a systematic way. 

In this work the system we are interested in consists of 
the adatom, the graphene, and the tip on the top of the 
adatom, as shown in Fig. [TJ As the hrst step, we explore 
the intrinsic equilibrium Kondo effect of the adatom as 
the conduction sea of electrons is limited to the graphene. 
The Kondo physics of the adatom generally depends 
on its explicit positions on the surface of the graphene. 
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When the adatom is absorbed on the top of carbon atom, 
the situation is quite simple since the orbital degree of 
freedom is not involved in the case. However, when the 
adatom is located at the center of the honeycomb, a high 
symmetry position, the adatom hybridizes with the car- 
bon atoms in both sublattices , 1 1 > 18 ~ — thus the orbital de- 
gree of freedom is involved which complicates the Kondo 
physics. For simplify, here we only focus on the first case 
and leave the latter one for the future study. When the 
Fermi energy is near the Dirac points, the Kondo effect is 
absent since a critical coupling strength must be satisfied 
for developing Kondo effect^ ' 16 ' 23 Once the Fermi energy 
is tuned to be away from the Dirac points, despite it is 
below or above the Dirac points, the Kondo resonance is 
found to exist in a wide parameter regime ranged from 
the Kondo regime to the mixed valence regime, even to 
the empty orbital regime. This is consistent with our in- 
tuition since the localized magnetic moment exists in a 
very wide parameter regime. When the system is in the 
Kondo regime, the Kondo resonance manifests as a sharp 
peak. This peak appears weaker and weaker when the 
impurity level becomes deeper and deeper. This is un- 
derstandable since the broadening of the impurity level 
becomes weak due to the low density of states of the 
graphene near the Dirac points. When the system is in 
the mixed valence or empty orbital regimes, the Kondo 
resonance develops, however, it does not show a sharp 
peak but a peak-dip structure, even an anti-resonance, 
which is a typical Fano resonance behavior . 24 ' 25 

In the next step, we consider a normal tip on the top of 
the adatom and the aim is to study the scanning tunnel- 
ing microscopy (STM) spectra of the Kondo resonance 
of the adatom in graphene. We first turn off the direct 
tunneling between the tip and the graphene. The Kondo 
resonance shows asymmetric behaviors that depends on 
the position of the Fermi energy, namely, it is above or 
below the Dirac points. Moreover, the Kondo tempera- 
ture strongly but asymmetrically depends on the explicit 
position of the Fermi energy which is in agreement with 
reported in the literature ! 23 ' 26 Subsequently, we turn on 
the direct tunneling between the tip and the graphene to 
study its influence on the lineshape of the Kondo reso- 
nance, which also is more realistic case. It is somehow 
surprising that the lineshape is almost unaffect by the 
additional direct channel. This result in fact originates 
from the unusual large asymmetry factor in graphene. 

The paper is organized as follows. In Sec. [TT] we 
take the single-impurity Anderson model to describe the 
adatom on the surface of the graphene and derive trans- 
port formula by Keldysh nonequilibrium Green's function 
theory. In Sec. IHII we solve numerically the equations ob- 
tained in a self-consistent way and discuss explicitly the 
Kondo effect in various situations including the cases of 
(i) the adatom plus the graphene, (ii) the adatom plus 
the graphene plus the tip, and (iii) the STM spectra with 
the direct channel between the tip and the graphene. Fi- 
nally, a brief summary is devoted to Sec. IIVI 
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FIG. 1. Schematic diagram of the STM measurement of a 
magnetic adatom absorbed on the top of a carbon atom site. 

II. MODEL AND FORMALISM 

The Hamiltonian of a magnetic adatom on graphene 
with an STM tip consists of the following terms 

H = H d + H g + H dg + Ht + H dt + Htg, (1) 

where 

H d = J2 £ ddU« + Udid a dU s , (2) 

describes the magnetic adatom with the impurity level e d 
and the on-site interaction U.— Here d\ (d a ) is a creation 
(annihilation) operator of localized electrons on magnetic 
adatom with spin er(=y, \) and a = —a. The second 
term in Eq. |T]) is the tight-banding Hamiltonian de- 
scribing graphene, which can be written as in momentum 
space 

H 9 = -tJ2(^al a b k(T + H.c?j, (3) 

ka 

a \ a ( a ka) and b\ a (bk a ) are creation (annihilation) oper- 
ators of electrons in sublattice A and B, respectively, t is 
the hopping energy between the nearest-neighbor carbon 

atoms, rjk = 2~^ =1 e lfe ' Fi with f\ = ax, r 2 = -f x+ ^^y, 
r 3 = — §a? — ^^y, and a « 1.42 A is the carbon- 
carbon spacing^ The low-energy transport properties of 
graphene are mostly determined by the nature of spec- 
trum around Dirac points K± at the corners of Bril- 
louin zone, where the dispersion is Ek ~ ±t>ir|fc|, the sign 
+ (— ) corresponds to the conduction (valence) band, and 
Vf = 3ta/2 (« 10 6 m/s) is the Fermi velocity of Dirac 
electrons. The tip Hamiltonian is described by a normal 
metal 

Ht = J2^tLtk„ (4) 

fcer 
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where t\ a (tka) represents the creation (annihilation) of 
an electron with energy on tip. The Hamiltonian 
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Vkbab 
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H.c. 



ka 



and 



H dt = E fe<4<A + H - c - 



(5) 



(6) 



ka 



represent electrons tunneling from magnetic adatom to 
graphene and to STM tip, respectively. Vkaa {Vkba) is the 
electron tunneling amplitude between magnetic adatom 
and sublattice A (B), and Vkta is the tunneling amplitude 
between magnetic adatom and STM tip. If the adatom is 
adsorbed above the carbon atom in the sublattice A, thus 
Vkaa = Vkga and Vkba = or in the sublattice B, thus 
Vkaa = and Vkba — Vkga- The direct tunneling Hamil- 
tonian between graphene and the STM tip is described 
by 



tg 



E t*. (*L« 

kk' cr 



H.c. 



(7) 



where Tt g is the hopping amplitude between the STM tip 
and graphene sheet. 



In steady state, the current from the STM tip to the 
graphene sheet can be calculated by time-dependent oc- 
cupancy number 



r= - 2 {{n, 

ka 



No 



N g )}) 



VkgaG^gka _ Vk*to-^a,tk<T (*»*', 



2e 



Re E r **s G 



< 

tk<7,gk' a 



(8) 



where JVj = E feCT *L^ and ^9 = E fcCT c L c ^- 
Here we use c\ a {cka) to denote the electron cre- 
ation (annihilation) operator in the graphene sheet. 
The current is related to the lesser Green's func- 



= *(*L(*'K(*) 



and Gf M , a {t,t') 



tions defined as G< tk(J {t,t') 

G<gka(t,t') = i (cl(t')da(t)) 

* ( c fc'cr (0 ^ These Green's functions can be calcu- 
lated by the equations of motion (EOM) approach based 
on nonequilibrium Green's functions theory. For a detail 
one can refer to Ref. 27f. For simplify, we assume that 
the tunneling amplitudes Vkga = V g and Vkta = Vt are 
independent on energy and spin. After some straightfor- 
ward derivations, one can obtain the current formula 



T = t E / ( £ ) - G ° t r * ( £ ) /* ( £ ) - r 9 ( £ ) u ( £ )i + ( r * ( £ ) - r « ( £ )) G ^ ( £ )} d£ 



+| / (/,(--) ■/.,(--)) 



n(e)-$( £ )£lm(G; ( £ )) 



de, 



(9) 



where G< r ' a (e) is the Fourier transformation of 
time-order Green's function G<' r ' a (t,t') defined 
by G<(t,f) - i (4(f) da (t)), G r a (a) (t,t>) = 
T^(±t T t')({da(t),di(t')}). T g (e) = n\V g \ 2 pg(e) 

and r t (e) = ir\V t \ pt{e) are the couplings between 
the magnetic adatom to the grarjhene and to the 
tip, respectively. n (e) = 7r|T tff | p t (e)p g (e) corre- 
sponds to hopping from the tip to graphene and vice 
versa. $ (e) = y/nU (e) T s (e) T t (e) is related to the 
interference between the different tunneling channels. 
p g (e) and pt (e) are the density of states of graphene 
and STM tip, respectively. The density of states for 
tip is simplified with a constant pt(s) = 1/21?, and 
the density of states for graphene can be evaluated 
with p g {e) = e(D - \e\)\e\/D 2 2L Here D is the half- 
width of conduction band. For convenience, we take 
r 9 (e) - r |e|/D where T = 7r\V g \ 2 /D. f g (e) and f t (e) 
are the Fermi distribution functions of graphene and 
the STM tip. The nonequilibrium Kondo effect can be 
discussed by differential conductance near the Fermi 



level. Thus, our final task to calculate the current is to 
solve the Green's functions involved. 

The equations for the retarded Green's functions and 
the lesser Green's functions can be derived by the 
Schwinger-Keldysh perturbation formalism^— which 
are given by, respectively, 

G r a (e)=g r d (s)(l + (([da,H I };dt)) r e ) , (10) 

G< {e)=gf{e) (l + ((K,Fz];4»3 

+g r d (e)(([da,H I ];d .))<, (11) 

where g r d {e) — 1/ (e — Ed + i0 + ) and g^ (e) — iTrf(e)S(e— 
Ed) are the noninteracting retarded and lesser Green's 
functions for the isolated magnetic adatom. In the 
above expressions the Hi = XJ d\d a d\d s + Hd g + H d t in- 
cludes the on-site Coulomb repulsion term on the adatom 
and the tunneling terms. While the first term will 
lead to higher-order Green's functions contained into 

(([do-, Hi] ;dt)}^ 0,< by using Zubarev notations^ the 
remaining two terms only lead to the hybridization effect. 
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Furthermore, the derivation of the EOM of the higher- 
order Green's function will also introduce more higher- 
order Green's functions, which have to be truncated in 
order to close the hierarchy of the EOM of the succes- 
sive Green's functions. Here we truncate this hierarchy 
of EOM by using Lacroix's scheme,— which is shown to 
be enough to capture the Kondo resonance, even at zero 



temperaturei 33 i 34 Since the derivation is standard we di- 
rectly write down the Green's functions in the infinite-U 
limit as follows 



Gl (e) 



1 - (n s ) - A s (g) 



g - e d - £ (e) (1 - A s (g)) - B s (g) ' 



(12) 



and 



G< (e) = G£(e)- 



(g) (1 - A, (g)) + B< (e) - A< (g) (G% (g)" 1 + £ (g)) 
1 - {n s ) - A a (e) 



(13) 



Here the self-energy S (e) = "^^rfc)*^"* » r ( £ ) = 
r s (g) + r t (g) and £< (e) = l %^ e \ e) , the average of 
occupation {n s ) = — i J ^^y-lm(G^ (g))efe, and 

, f v i y g( £ , £ Q gi ( £ o+zn( £ Q)(Gt( £ or 

ttJ l + 7rn(e') 

(14) 



7T J 1 + 7rn(e') v ; 



A< _ gi (e) Re (G| ( £ )) + i» (g) Im {G r s (g)) ^ 



1 + Trn (g) 



B< (g) = iSi (g) + i (Sa (g) + S 3 (g)) Im (G^ (g)) 



$(g)Im(G; OQ) 
1 + Trn (g) 



(17) 



where g (g, g') = 1/ (g - g' + i0+), fi (g) = $ (g) / 9 (e) + 
($ (g))* /t (e), S^g) = / 9 (g)r g (g) + f t (e)T t (e), 

H 2 (g) = / ff (g)r s (g)r 9 (g) + / t (g)r t (g)r t (g), 
Ss(e) = r 9 (g)r t (g)(/ 9i (g) + / t9 (g)) with 



fgt(e) = / 9 (g)(l-/ t (g)). 



At the hand of the 



above expressions, we can calculate numerically the 
Green's function G r a (g) and (g) = (G r a (g))* in a 
self-consistent way, and further the lesser Green function 
G^ (g). Finally, we obtain the current I by Eq. ©. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In the following numerical calculation, several input 
parameters are tuned in order to study the system in dif- 
ferent physical situations. The atomic level ed is changed 
to study the different regimes ranged from the Kondo 



regime, to the mixed valence regime, even to the empty 
orbital regime^ The chemical potential of graphene 
can be tuned by gate-voltage, which makes the Fermi 
level is located at (/i = 0), above (fi > 0) or below 
(fi < 0) the Dirac points. The different chemical poten- 
tial will lead to quite different physics, as shown below. 
The voltage bias V s d is also applied between STM tip and 
graphene sheet to calculate the differential conductance. 
For convenience, we take To as the units of energy in the 
calculations. As a systematic study, we first to calculate 
the Kondo effect of the adatom plus the graphene when 
the STM tip is omitted. The aim is to explore the in- 
fluence of the linear dispersion around the Dirac points 
to the Kondo effect and compare it with the adatom on 
normal metal surface. Subsequently, we consider the in- 
fluence of the STM tip and explore the Kondo resonance 
when both the normal metal and the graphene are cou- 
pled to the adatom. Finally, we turn on the direct chan- 
nel between the STM tip and the Graphene and check its 
effect in the tunneling spectroscopy, which is more closely 
related to the experiment. 



A. Kondo effect of the adatom without the STM 

tip 

In order to explore the effect of the linear dispersion 
around the Dirac points, we first the system composed 
of the magnetic adatom plus the graphene and omit the 
STM tip. As pointed out by Uchoa et al.,- 3 it is more 
easier to form the local moment for the magnetic impu- 
rity on the surface of the graphene than in the normal 
metal host. This conclusion has a direct consequence 
that one can observe the Kondo resonance in a more 
wider parameter region for the former than the latter 
one. This is really true as shown in Figs. [5] and [3] where 
we show the local density of states of the adatom given by 
Pd( £ ) — — ^ImG^e). In Fig. [5] we consider \i = 0.5r , 
namely, the Fermi level is above the Dirac points. We 
change the atomic level from the empty orbital regime, 
as shown in Fig. [^a) to the Kondo regime [see Fig. [2ff)]> 
around the Fermi level the local density of states show 
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FIG. 2. The Kondo resonance of the magnetic adatom on 
graphene with different atomic levels ranged from the empty 
orbital regime to the Kondo regime. From (a) to (f) the 
atomic level is e d = 2.5T , 1.5F , 0.5F , 0, -0.5r , -1.5r . 
The inset in each panel is enlarged view around the Kondo 
resonance. The chemical potential is /i = 0.5Fo, the temper- 
ature T = 0.00005r , and the half-bandwidth D = 5r . 




FIG. 3. The same as Fig. |2]but for fx = -0.5r . From (a) to 
(f) the atomic level is e d = 1.5r , 0.5r , 0, -0.5r , -1.5r , 
~2.5r . 



special feature, i.e., the Kondo resonance, although their 
shapes are quite different. In the Kondo regimes, the 
Kondo resonance exhibits as a sharp peak around the 
Fermi level but one should notice that with increasing 
the atom level the peak becomes gradually asymmetric, 
which is quite obvious in Fig. WLA) for Ed = 0. Continuing 
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FIG. 4. Comparison between the density of states around 
the Fermi level in the empty orbital regime for the metallic 
(red dashed lines) and the graphene (solid lines) hosts, (a) 
fx = 0.5r , E d = 2.5r , and (b) fx = -0.5T , e d = 1.5r . 



to change the atomic level to Ed = 0.5Fo, the Kondo reso- 
nance has a typical peak-dip structure. Finally, we move 
the atomic level to the empty orbital regimes, the Kondo 
effect manifests as an anti-resonance, which is clearly seen 
from the inset of Fig. E^a) and (b). These features are 
nothing but the Kondo resonance. They can be under- 
stood based on our previous worki^ In such a system, 
the Kondo resonance results from the coherent supposi- 
tion between the Kondo peak and the broadening of the 
impurity level due to the hybridization with the conduc- 
tion electron. In the Kondo regime, the broadening of 
the impurity level is quite weak around the Fermi level 
so the Kondo resonance shows as a sharp peak. When 
the system approaches to the mixed valence regime, the 
broadening of the impurity level becomes more and more 
sizeable and as a result, the Kondo peak distorts grad- 
ually due to the interference, and finally it becomes an 
anti-resonance when the system is in the empty orbital 
regime since in this case the Fermi level is located at the 
left-hand side of the broadening impurity level. The es- 
sential reason that the Kondo effect can exist in a wide 
parameter region is the impurity level is broadened exten- 
sively due to the linear density of states of the graphene 
around the Dirac points. Therefore, the adatom on the 
surface of the graphene provides an idea system to study 
the Kondo effect in various parameter spaces. The same 
physics has been observed when the Fermi level is below 
the Dirac points, as shown in Fig. |3][a)-(f). 

To confirm above observation, we compare the local 
density of states around the Fermi level for the adatom 
on graphene and on metal surface in the empty orbital 
regime, as show in Fig. 01 For metal host, one cannot 
observe any feature, which is in contrast to that in the 
graphene, where an obvious anti-resonance exists. The 
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FIG. 5. The Kondo resonance around the zero bias varies 
with the chemical potential (a) fj, < and (b) fi > 0. (c) 
The Kondo temperature Tk versus the chemical potential p. 
The other parameters used are = — ro where Ft = O.lFo, 
T tg = and D — 5Fo. (d) shows schematically the cotunnel- 
ing processes with the chemical potential p < and n > 0, 
respectively. 
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without the direct tunneling T tg channel for two cases (a) 
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Ttg = 0.2ro. The other parameters used are the same as Fig. 
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results are true either the chemical potential is above or 
below the Dirac points. 

B. Kondo Effect of the adatom with the STM tip 

When the STM tip is considered, the adatom couples 
with two baths, one is the STM tip which is a normal 
metal with a constant density of states, as studied in the 
conventional Kondo physics. The other is the graphene 
which has a linear density of states around the Dirac 
points. We first study the nonequilibrium Kondo effect 
in this case. Therefore, we omit the direct channel be- 
tween the STM tip and the graphene. In the calculations, 
the coupling between the adatom and tip is taken with 
T t = 0.ir . Fig. [5] shows the differential conductance 
around the zero bias. When the chemical potential /i is 
tuned to meet the Dirac points, the curve at the zero 
bias shows a kink but no visible Kondo peak, which is 
a consequence of the zero density of states at the Dirac 
points. When the chemical potential fi is tuned to be 
away from the Dirac points, either above or below the 
Dirac points, the zero bias peak occurs, which shows the 
Kondo effect. When the chemical potential below the 
Dirac points, the more larger is, the more wider the 
peak is. This is easy to understand since the larger [i 
means the larger density of states at the side of graphene. 
When fi is above the Dirac points, the width of the Kondo 
resonance keeps almost unchange (see Fig. HJc)), which 
also can be obtained from the chemical potential depen- 
dent Kondo temperature^ 2 - However, when /i is below the 
Dirac points, the Kondo temperature Tk shows obvious 
dependence on [x. As a result, the Kondo temperature 



is obvious asymmetric, which is consistent with that re- 
ported in the literature . 23 ' 26 This is due to the fact of 
the particle-hole asymmetry in graphene. In addition, 
one also notes that in both cases the curves show a kink 
when the bias scans through the Dirac points. 

C. Scanning tunneling spectroscopy with direct 
channel 

In this subsection we discuss the scanning tunneling 
spectroscopy of the adatom on the surface of graphene 
by switching on the direct tunneling channel between the 
STM tip and graphene, a more realistic case in experi- 
ments. Due to the additional direct channel, one may ex- 
pect that the STM spectroscopy should be changed dra- 
matically leaded by the interference between the direct 
and the indirect channels, as usually observed in a system 
with the adatom on the surface of the normal metals.— 
However, the result is somehow surprising, as shown in 
Fig. [6] In comparison to the spectroscopy without the 
direct channel, except for the amplitudes, the shape al- 
most keeps unchanged for both fi > and \i < 0. To 
check the underlying physics of this result, let us recall 
the Fano resonance theory^ In the STM case, due to the 
adatom scattering and its interference with the direction 
channel, the differential conductance of the conduction 
electron can be written as2£ 

^(x(q 2 - l)Im G(e) - 2qRe G(e), (18) 

where G{e) denotes the adatom Green's function and 
q ?s — Re g(e)/lmg(e) is the Fano asymmetry factor. 
Here g(e) is the graphene Green's function with Re g(e) = 
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(e/D 2 )\n\e 2 /(D 2 - e 2 )\ and lmg{e) = -n\e\/D 2 Q{D - 
\e\). Thus the q factor is 2 sgn(/x)/7r In \fi\/D near the 
Fermi levels For graphene, D ~ 6eV (Ref. [36]) and 
if \fi\ maximally is in the order of O.leV, so |g| > 6. 
That means that the first term in Eq. (fT5|) dominates 
the second term, as a result, the dl/dV has the almost 
same profile as ImG(e), the local density of states of the 
adatom. This is the reason that the shapes of the Kondo 
resonance around the zero bias is almost same despite 
the direct channel is switched on or off. 



IV. SUMMARY 

In summary, we study systematically the Kondo effect 
of a magnetic adatom on the surface of graphene and its 
STM spectra. The main conclusion is that the Kondo res- 
onance can exist in a much wide parameter region ranged 
from the Kondo regime, to the mixed valence regime, 
even to the empty orbital regime. In the mixed valence 
or the empty orbital regimes the Kondo resonance show 
specially lineshape due to the Fano resonance since in 



these regimes the broadening of the impurity level pro- 
vides a significant background. The STM spectra have a 
particle-hole asymmetry when the chemical potential in 
graphene is tuned by the gate voltages. Furthermore, it 
is found that the direct tunneling channel between the 
STM tip and graphene has no obvious effect on the line- 
shape of the Kondo resonance around the zero bias, which 
is a consequence of the unusual large asymmetry factor 
in graphene host. The rich behaviors of the Kondo res- 
onance obtained indicate that the graphene is an ideal 
platform to study the Kondo physics and these results 
are useful to further stimulate the relevant experimental 
works in this system. 
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